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Numerical simulations of Ramsey fringe formation in double Bose condensates are carried out in 
two spatial dimensions. The effects of meanfield nonlinearity and diffusion in the condensates give 
rise to new features in the Ramsey fringes, namely spatial variation and an asymmetry between 
positive and negative field detunings. We introduce the concept of an effective detuning, which for 
moderate interpulse times provides a qualitative understanding of the new features. 



I. INTRODUCTION 

Recent experiments carried out by the JILA and MIT 
groups U' o have shown the rich behaviour that can 
occur in double condensate systems. In particular the 
JILA group have used Ramsey's 61 separated pulse tech- 
nique to produce interference fringes which demonstrate 
the long term phase coherence of their condensates. Al- 
though Ramsey's original theoretical description explains 
qualitatively the behaviour that is seen, additional mech- 
anisms present in condensates, namely nonlinearity from 
the mean field collisions and diffusion due to kinetic en- 
ergy, complicate the condensate case. In this paper, we 
present a theoretical study of Ramsey fringe formation in 
two spatial dimensions, including the effect of the mean 
field collisions and kinetic energy. The JILA experiment 
minimises the nonlinear effects by only partially overlap- 
ping the condensates, so that the experiment is carried 
out in the low density (linear) region. Here, we consider 
the case where the condensates are fully overlapped, and 
show that the additional condensate mechanisms can give 
rise to spatial variation of the fringes across the conden- 
sate, and to an asymmetry between positive and neg- 
ative detunings. Eventually, for large separation times 
between pulses, the regularity of fringes predicted by 
the original theory is disrupted by the additional mech- 
anisms. 



II. FORMULATION 

In Ramsey's original work @ a beam of atoms in the 
ground state was passed through two separated oscilla- 
tory field regions. The interaction with the field in the 
first region transfers some population into the excited 
state and creates a magnetic dipole. Between the oscil- 
lating fields the atomic populations remain fixed, but the 
dipole undergoes free precession. The interaction with 
the second field causes an additional transfer of popula- 
tion, which is sensitively dependent on the total phase 
accumulated by the dipole in the field free region. If the 



pulses are separated by a fixed time T, the measured pop- 
ulation shows a fringe pattern (as a function of the oscil- 
lating field frequency) known as a Ramsey fringe, which 
is analogous to that obtained in a 2-slit interference ex- 
periment. The experiment can also be carried out with a 
fixed field frequency and varying pulse separation times 
T, which leads to a similar interference pattern, with the 
population varying according to cos(AT), where A is the 
atom-field detuning. 

In the double BEC system, the two atomic states of 
the original Ramsey system become the two conden- 
sate states of different internal states, and the oscillating 
fields are 2-photon ir/2 pulses (of duration r = 7r/(2fi), 
where ft = + A 2 , and fin is the 2-photon Rabi 

frequency u ) . We choose in this paper to scan the pulse 
separation time T and hold the detuning constant—fol- 
lowing the procedure used in the JILA experiment 

In the meanfield limit, the double condensate system 
can be described by the coupled Gross-Pitaevskii (GP) 
equations o< O, 
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where the spatial and temporal coordinates, r and t , are 
scaled as in Ballagh et al 0. The nonlinearity param- 
eter C is proportional to the total number of atoms N 
in the condensate and the intraspecies s-wave scattering 
length, while the factor k indicates the relative trapping 
strength of state |2). In this paper, we take k = 1, so that 
both components are equally trapped, as in the JILA 
experiment. The factor w gives the ratio of inter- to 
intra- species scattering length, and includes the effect 
of wavefunction symmetry. We have assumed that the 
intra-species s-wave scattering length an is the same for 
each internal state of the atom, so that an = 022- This is 
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reasonable for the alkali atoms, and makes the C values 
the same in both (|) and (§), but our results could eas- 
ily be generalised to situations of unequal intra-species 
scattering length. Growth and loss processes are ignored 
here, so that N is conserved. The coherent coupling be- 
tween condensate components is described by the terms 
containing V and A , where V is half the bare Rabi fre- 
quency, and A is the detuning of the coupling field from 
the bare atomic resonance 0. 

It is useful at this point to introduce a new quantity 
A e f f , which is the difference between the diagonal terms 
in Eqns (0) and @, 

A e// (r,i) = C(w - l)(|ViM)| 2 - |^ 2 (r, t)\ 2 ) + A. (3) 

This quantity, which we will interpret as an effective de- 
tuning, will prove helpful in interpreting the effect of 
the nonlinear mechanism on fringe formation. It is clear 
that when w = 1 (all scattering lengths the same), or 
\tpi(r,t)\ 2 w | -02 (r, i) | 2 (near the edges of the wavefunc- 
tions), A e ff(r, t) reduces to the usual value of the detun- 
ing A. 



III. RESULTS 

Equations (Q) and ^ are solved for a succession of T 
values at fixed detuning A, using a modified split-step 
Fast Fourier Transform method. Initially all the pop- 
ulation is in component 1, which we choose to be an 
eigenstate of the uncoupled single component system. A 
grid of 512x512 points is used for the simulations, over a 
spatial range of —20 < x, y < 20. In our results, we have 
chosen the case of w — 2, which makes apparent the role 
of the meanfield nonlinearity, and is realistic for sodium. 

A typical result is shown in Fig. 0(a), where the central 
condensate density (i.e. at x = 0, y = ) of each com- 
ponent immediately following the interaction with the 
second oscillatory field is plotted as a function of the 
interpulse time T. The figure displays the characteris- 
tic Ramsey fringe pattern, however more careful inspec- 
tion of our results shows that the additional condensate 
mechanisms cause significant modifications to the fringe 
behaviour. This is most obviously apparent in the spatial 
density distributions of each component, and in Fig. ^ we 
plot the full two dimensional behaviour for the same case 
as Fig. 0(a), at the interpulse time of T = 0.6674. Spa- 
tial structure is now very evident in component 2, and 
this can not be accounted for in the simple Ramsey the- 
ory. Density slices taken across the x axis allow a more 
quantitative evaluation of the spatial structure, and in 
Fig." 
Fig. 



2 we plot, for the same simulation parameters as in 
1(a), a sequence of density slices corresponding to a 
sequence of equally spaced interpulse times. It is clear 
from this sequence that component 2 does not replicate 
the initial eigenstate of component 1, but instead devel- 
ops a dip in its centre, which increases as T increases. 



A further significant difference of the condensate case is 
seen when we consider the effect of changing the sign 
of the detuning A. Ramsey's original theoretical expres- 
sion for the fringes is exactly symmetric to the change 
A — > —A. In the condensate case, when we repeat the 
previous calculations with the opposite sign for detuning 
(and all other parameters unchanged), we obtain signifi- 
cantly different results, as shown in Fig. |l|(b), and even 
more dramatically in Fig. ^. 

Comparing Figs. |l|(a) and (b). we see that the period 
of the Ramsey fringe at the centre of the wavefunction 
is distinctly different for detuning of opposite signs: in 
Fig. |l|(a) (where A = 10) the first minimum of compo- 
nent 1 occurs for T = 0.47, while in Fig. @(b)) (where 
A = —10) the first minimum of component 1 occurs for 
T = 0.61. Comparing Figs. || and f|, for which the se- 
quence of times are identical and only the sign of the 
detuning has changed, the contrast is immediately ap- 
parent. In Fig. |[ component 2, which is decreasing as 
a whole during the sequence, develops a central dip. In 
Fig. |i| on the other hand, component 1 is increasing as a 
whole, and develops a central dip. 

These results can be understood qualitatively in terms 
of the effective detuning A e f / (r, t) , and the well known 
Thomas Fermi approximation. The effective detuning 
A e ff(r,t) is a dynamic quantity that depends on the 
instantaneous population density difference \tpi(r,t)\ 2 — 
1^2 (r, t)\ 2 . Within the Thomas Fermi approximation 
(which neglects the V 2 terms in the GP equation) no 
mixing occurs between different spatial points, and cal- 
culations are greatly simplified. For example the atomic 
dipole created by the first pulse can be calculated, with 
the full dynamic behaviour of A e f f (r, t) taken into ac- 
count, by solving an ordinary differential equation. How- 
ever, good estimates of the behaviour of the dipole can 
be obtained even more simply, as we will now discuss. 
We note first that during the interval between pulses, the 
populations at position r do not change (as a result of ne- 
glecting the V 2 term), and the phase accumulated by the 
atomic dipole is determined by the value of A e ff(r, r) at 
the end of the first pulse (which is of length r) . The value 
A e jj(r, t) is determined by the inversion density at time 
r, and for a constant detuning A the inversion following 
the first pulse is given by the well known formula o , 
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where wq is the initial inversion and Qq = V 'AV 2 + A 2 . 
At the edges of the condensate, where the densities 
are low, A e // is essentially A, and thus the first ir/2 
pulse (OqT = 7r/2) produces a resulting inversion of 
wo(A/flo) 2 . Elsewhere across the wavefunction, we can 
approximate the effect of the dynamically changing de- 
tuning by using a temporal average of A e //(r, f) during 
the pulse, which we will write as A e //(r). At the centre 
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FIG. 1. Condensate density of component 1 (solid line) and 
component 2 (dashed line) at the trap centre, as a function 
of interpulse times T. (a) A = 10 (b) A = -10. The density 
is taken immediately following the second pulse. Parameters 
are V = 10, G = 200, w = 2, r = 0.0703 



(a) 




FIG. 2. Spatial distribution of condensate density of (a) 
component 1 (b) component 2, for an interpulse time of 
T = 0.6674. Parameters are as in Fig. FH(a) 
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FIG. 3. Condensate density of component 1 (solid line) 
and component 2 (dashed line) along the x axis, for a se- 
quence of interpulse times (a) T = 0.5971 (b) T = 0.6322 (c) 
T = 0.6674 (d) T = 0.7025. The detuning is A = 10, and all 
other parameters are as in Fig. Note that (c) corresponds 
to Fig. | 
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FIG. 4. Condensate density of component 1 (solid line) and 
component 2 (dashed line) along the x axis, for a sequence of 
interpulse times (a) T = 0.5971 (b) T = 0.6322 (c) T = 0.6674 
(d) T = 0.7025. The detuning is A = -10, and all other 
parameters are as in Fig. |IJ 



0.03 




12 3 

interpulse time 

FIG. 5. Condensate density of component 1 (solid line) and 
component 2 (dashed line) at the trap centre, as a function 
of interpulse times T. All parameters are as in Fig. |l|(a) 

of the condensate, the densities are largest, and the aver- 
age over time of |?/>i(r, t)\ 2 — \4>2(r, t)\ 2 is positive during 
the pulse, so that |A e //(r)| exceeds |A| (for positive A). 
Thus from Eq.([|) we see (in the regime where V < A) 
that the change in the initial inversion density is less at 
the centre than at the edges. Thus for positive detuning, 
A e f f (r, r) is greater at the centre of the wavefunction 
than at the edges, and thus the subsequent phase accu- 
mulation between pulses is greater in the centre. We thus 
expect that the Ramsey fringes will have a greater fre- 
quency in the centre of the wavefunction, which provides 
the explanation for the behaviour in Fig. [|. There we 
see that population transfer out of component 2 is more 
rapid at the centre than at the edges, which means the 
interpulse minimum is reached first at the centre. In the 
case of negative A, the magnitude of the effective detun- 
ing will be less at the centre than at the edges [provided 
| A | is larger than | C (w ~ 1 ) ( | fa (r , t) \ 2 - | Xb 2 (r , t) | 2 ) | ] . This 
means the frequency of the Ramsey fringes will be less 
at the centre of the wavefunction than the edges, and ex- 
plains why, in Fig. [|, the edges of component 1 reach their 
maximum more rapidly than the centre. It is also now 
clear that Ramsey fringes at the centre of the wavefunc- 
tion will have a greater frequency for positive detuning 
than for negative detuning, which indeed is the behaviour 
that is revealed by comparison of Figs. 0(a) and (b). 

Finally we present, in Fig. ||, the behaviour of the Ram- 
sey fringes for a large range of interpulse times T.At long 
times, the effects of diffusion combined with meanfield 
nonlinearity produce irregular periodicity. We note too 
that a plot of the spatial dependence of the population 



densities at a given large value of T will usually show sig- 
nificant modulation. Once diffusive mixing becomes sig- 
nificant, the predictions of our effective detuning, which 
relies on the Thomas Fermi approximation, are no longer 
valid. The fringe behaviour must then be obtained from 
a full solution of the coupled GP equations (|l|) and (||). 



IV. CONCLUSIONS 

We have carried out numerical simulations in two spa- 
tial dimensions of Ramsey fringe formation in double 
Bose condensates. Our calculations include the effects 
of both nonlinearity and diffusion, and in contrast to the 
reported JILA experiment, we have considered the case of 
fully overlapping condensates, which enhances the effect 
of nonlinearity. The mechanisms of diffusion and non- 
linearity give rise to new features in the Ramsey fringes, 
namely spatial variation of the fringes, and an asymme- 
try between positive and negative field detunings. By 
introducing the concept of an effective detuning we have 
obtained, for moderate interpulse times, a qualitative un- 
derstanding of the new features. However the effective 
detuning accounts only for nonlinear effects and eventu- 
ally, at long enough interpulse times, diffusion becomes 
important and the full GP equations must then be solved. 
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